Test of a LYSO matrix with an electron beam between 100 and 500 MeV for KLOE-2 
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Abstract 

The angular coverage extension of the KLOE-2 electromagnetic calorimeter, from a polar angle of 20° down to 8°, will increase the 
C^) multiphoton detection capability of the experiment enhancing the search reach for rare kaon, 77 and 77' prompt decay channels. The 
basic layout of the calorimeter extension consists of two small barrels of LYSO crystals readout with APD photosensors aiming to 
£H achieve a timing resolution between 300 and 500 ps for 20 MeV photons. The first test of a (5.5x6x13) cm 3 prototype for such a 
^ detector was carried out in april 2009 at the Beam Test Facility of Laboratori Nazionali di Frascati of INFN with an electron beam 
'™~ 5 from 100 to 500 MeV. In the selected energy range, we measured a light yield of 500H-800 p.e./Mev, an energy resolution which can 
be parametrized as 0.05 ©0.01 /(E/GeV)© 0.015/ VE/GeV, a position resolution of 2.8 mm and a timing resolution of 200n-300ps. 
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a 

. ^ Crab-waist and a large Piwinsky angle has been proposed and 
^ tested U] to improve the reachable luminosity at the Frascati 
J"^* ^-factory, DA<DNE, an e + e~ collider running at the center of 
Q^mass energy of the <\> resonance. The success of this test mo- 
'— 'tivated the startup of a new experiment, named KLOE-2 [|2|], 
j which aims to complete the KLOE Q3D physics program and 
perform a new set of interesting measurements. The first run- 
■^j- ning phase will start at the end of 2009 with the goal of collect- 
ing ~ 5 fbr 1 in one year. For this first run, the modification of 
1— ' the KLOE detector will be minimal with the only introduction 
' \ of a 77 tagging system. A second phase, for a longer data tak- 
ing, will require another set of upgrades all concentrated around 



1. Introduction 

In the last years, a new machine scheme based on the 



the beam-pipe consisting of: (i) an inner tracker, IT, based on 



^ Cylindrical GEM technology, (ii) a tile calorimeter surrounding 
• • the inner quadrupoles , QCALT , and (iii) a calorimeter between 
. ^ e interaction point, IP, and the first inner quadrupole. In the 
following, we describe the work done to design and test the 
$— 1 possibile solution for the latter detector. 
. 5^ . In Fig. Q] we show a zoomed-view of the available region 
around the IP which can be used to extend the angular cover- 
age of the main electromagnetic calorimeter, today limited at 
a polar angle of ~ 20°, with the addition of a new dedicated 
calorimeter. Assuming to be able to lower the minimum polar 
angle for photon detection down to 8°, this will enhance the 
multiphoton detection capability of the detector for the search 
of rare decays of kaons, 77 and 77' mesons. As an example, in 
2007 the KLOE experiment has measured a BR(^fs — > 77) J3l 
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Figure 1: Zoomed-view of the IP region. The area available for the new 
calorimeter lies between the inner sphere and the closest quadrupoles. 



in good agreement with 0(p 4 ) calculation of Chiral Perturbation 
Theory but which differs of three standard deviations from the 
same measurement carried out by the NA48 experiment. The 
proposed upgrade will allow to improve the S/B ratio of a fac- 
tor of three thus allowing to solve this puzzle. Other searches or 
branching ratio determinations, such for instance K s — > 3n° or 
77 — > 7r°77, will directly benefit from this extension. The only 
available area to place a calorimeter lies between the end of the 
spherical beam pipe, of 10 cm radius, and the first quadrupole, 
positioned at 30 cm distance from the IP. 

2. CCALT: a Crystal Calorimeter with Time 

The discussion of the previous section indicates that this 
calorimeter has to be very dense, with a small value of radia- 
tion length, Xo, and Moliere radius, R m , not hygroscopic and 
with a large light output to improve photon detection efficiency 
at low energy (from 20 to 500 MeV). Moreover, the calorime- 
ter has to be extreemely fast in order to allow prompt photon 
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Figure 2: Composition of crystals used for the matrix prototype. 



reconstruction in an environment with a large background rate 
(~ 1-r 5 MHz) of secondary showers generated by off-axis e + , 
e~ coming from intra-bunch scattering (Touschek effect). Pre- 
liminary simulation studies indicates the need to reach a time 
resolution of 300-H 500 ps for 20 MeV photons. 

A suitable solution is offered by a crystal calorimeter with 
good timing performances, named CCALT. A first detector lay- 
out consists of two concentrical barrels of 24 crystals each, with 
transversal dimensions of 2x2 cm 2 and longitudinal length be- 
tween 13 and 15 cm. The best crystal choice matching the re- 
quirements is provided by the new generation of Cerium doped 
Lutetium Yttrium Orthosilicate, LYSO, which has Xq and Rm 
values (1.1 and 2 cm) comparable to the ones (0.9 and 2 cm) of 
the Lead Tungstanate, PbW04, with the advantage of a much 
larger light yield (x 300). On the negative side, LYSO shows a 
scintillation emission time (tlyso — 40 ns) slower than PbW04 
(jpbwo = 10 ns). However, from the basic scaling law of the 
timing resolution, cr, = r/ y"Np.e., we expect LYSO to be a 
factor four more performant than PbW04. 

In the final location of the CCALT inside KLOE-2, the pres- 
ence of an axial magnetic field of 0.52 kGauss forces the usage 
of silicon based photodetectors. Due to the high photon yield, 
the readout with APDs is a valid solution since, at the lowest 
photon energy of 20 MeV, the collected photoelectrons will be 
around ~ 10.000 which corresponds to ~ 12 pC assuming an 
average gain of 300 and an amplification stage of x25, well 
matching the ADC sensitivity of the KLOE calorimeter (100 
fC/count). In the following, we specifically considered only 
the Hamamatsu S8664-55, which has an active area of 0.5x0.5 
cm 2 , fast timing characteristics and a quantum efficiency be- 
tween 65 and 85% in the wavelenght range of interest (390~ 
500 nm) for the LYSO emission spectra. 

In march 2009, we have built a medium size crystal matrix 
prototype with transversal radius larger than 2 R m , longitudinal 
dimensions being constrained by budget limits to be between 
13 and 15 cm (corresponding to 11 -f 12 Xq of longitudinal 
containment). The prototype consists of an inner matrix of 10 
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Figure 3: Distribution of Qtot for single electron events at 100 MeV (top) with 
a logG fit superimposed and at 500 MeV (bottom) with a gaussian fit. 



LYSO crystals readout by APD and an outer matrix, for leak- 
age recovery, composed by 8 PbW04 cristals readout by stan- 
dard Hamamatsu Bialcali photomultipliers of 1,1/8" diameter. 
To test the quality of the crystals offered by different vendors, 
the inner matrix has been assembled in three rows (see Fig. [2]) 
composed (from bottom to top) as follows: 

• 3 LFS crystals from Zecotek of 2x2x13 cm 3 , 

• 2 LYSO St.Gobain crystals of 1.5x1.5x15 cm 3 , 1 LYSO 
St.Gobain + 1 LYSO Seionix crystals of 1.5x1.5x13 cm 3 



1 LYSO St.Gobain crystal of 2x2x15 cm 3 , 
Seionix crystals of 2x2x13 cm 3 . 



2 LYSO 



The LFS from Zecotek is a Luthetium Fine Silicate crystal, 
with very similar properties to LYSO. Each crystal is wrapped 
with 100 jum of tyvek on the lateral faces, leaving free both 
the front and end faces, thus allowing to bring calibration light 
pulses through an external LED and a fast change of the pho- 
tosensors readout. Each APD is inserted in a PVC mask with 
the amplifier soldered on its anode and then mecchanically po- 
sitioned inside a stainless steel box closed by a PVC cap with 
only the electronic pins coming out for connection to HV and 
readout cables. An external holder takes the PMs in position 
for the readout of the outer crystals while allowing to press the 
boxes containing the APDs. The optical connection of the pho- 
tosensors with the crystals is done with optical grease. The 
amplifiers are based on the MAR8A+ chip from Minicircuits, 
with a gain factor of 25 and a bandwidth of 1 GHz. 
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Figure 4: Shower example generated by a single event of 100 MeV electron in 
the matrix as shown by a full simulation based on Geant-4. 

3. Test results with electron beams. 

We have taken data at the Beam Test Facility, BTF, of LNF 
for two weeks in april 2009. The matrix was positioned at the 
center of the beam axis with an area delimited by a cross of two 
finger BC408 scintillators of 1 x 0.5x5 cm 3 dimensions, /l,/2. 
In most of the tests, the fingers were aligned in such a way to 
define a beam spot of 1 x 1 cm 2 . In front of the fingers it was 
also present a beam position monitor, BPM, of the BTF group, 
consisting of sixteen horizontal and vertical scintillator strips 
readout by two Multi Anode PMs. Each strip is built by three 
1 mm diameter scintillating fibers thus providing an accuracy 
below 1 mm on the beam localization. 

We have triggered by using a replica of the spill signal from 
the Linac adjustable from remote in order to correctly put the 
signals in time. Prototype signals were splitted and discrimi- 
nated by means of the standard KLOE calorimeter SDS boards. 
We acquire data with the KLOE-2 daq system, VME based, 
reading out KLOE ADC and TDC boards with a sensitivity of 
100 fC/count and 50 ps/count respectively. 

The BTF beam allows to set the average number of electrons, 
< Ne > arriving to the detector by adjusting dedicated collima- 
tors. However, we could not lower < Ne > below 1 for all 
beam energies and we cleaned the data offline by cutting on 
pulse height on the fl, f2 counters. 

Observing the response of the prototype to single electrons, 
we realized that the outer matrix was not properly working due 
to an unexpected optical cross-talk between crystals. This was 
discovered by pulsing a UV LED (at 380 nm) in the front face 
of each LYSO crystal to let it scintillate and testing, both at the 
scope and with a fast data taking, the amount of signal observed 
on the contiguous crystals. We observe large cross-talk only on 
the PbW04 of the outer matrix. We believe this to be due by 
a cooperation between a light leak through the tyvek and the 
different amplification gains between PMs and APDs. In the 
following, we therefore report only results related to the inner 
matrix. 

By using the UV LED, we have first equalized each channel 
at 10% level by proper HV adjusting. We have then calibrated 
the calorimeter response of each channel with minimum ion- 
izing particles, m.i.p., crossing the calorimeter hortogonally to 
the crystal axis. To do this, we have collected cosmic ray runs 
with a dedicated trigger defined by two plastic scintillators po- 




Figure 5: Dependence of the energy resolution on beam momentum 

sitioned one above and one below the prototype. We get cr pe d 
of 5 counts and a m.i.p. peak, M,, of around 100 counts for the 
smaller size crystals. The statistical precision on the peak de- 
termination is ~ 1 %. The total response of the detector is then 
defined as: 

Qtot = J] (Qi - Pd x Mo/Mi, (1) 

where Qi and P, are the collected charge and the pedestal of the 
i-th channel, Mo represents an average calibration of all chan- 
nels in counts and the calibration for the larger crystals is cor- 
rected for the different size. In Fig. [3] we show the distribution 
of Qtot for a beam of 100 and 500 MeV respectively after having 
selected single electron events with a cut on the finger scintilla- 
tors. A surviving fraction of events with more than 1 electron is 
still observed in the matrix, expecially at low energies. We have 
fit the distribution corresponding to one electron either with a 
simple gaussian, centered around the peak, or with a logarith- 
mic gaussian, logG, function as follows: 

N ■ eX p(-JL/„(l _ JL {E - E peak )f - ^) (2) 

where N is a normalization factore, r\ represents the asymme- 
try, Epeak the most probable value of the distribution, ctq = 
Y^rsinh~ x (2. 36/7/2) and <r E = is the resolution. 

By performing a linear fit to the distribution of E pea k vs Pheam, 
we get a slope of 7.5+ 0. 1 counts/MeV which sets the M-, value 
to be ~ 16 MeV for a small crystal consistently with an expected 
energy loss of ~ 10 MeV/cm. At the running voltages of 410 V, 
the expected APD gain varies between 300-500 from which we 
estimate the light yield to be between 500 and 800 p.e./MeV. 

To understand the different terms of the energy dependence 
of the energy resolution, we are earring out a full simulation of 
the prototype based on Geant-4. In Fig. [4] we report one ex- 
ample of 100 MeV electron shower as seen from the MC pro- 
gram. We have also simulated the dead space of 100 /urn of 
tyvek between crystals and activated the optical transportation 
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Figure 6: Distribution of Qtot for events without electrons for runs with (up) 
500 MeV, (down) 100 MeV beam. 

of photons. Studies on photoelectron collection efficiency are 
still underways. This simulation indicates that there is a large 
leakage term between 5 and 4 % from 100 to 500 MeV. 

In Fig. [5] we show the energy dependence of the energy res- 
olution measured on data which has been fit with the following 
equation: 

cr E IE = a® b/(E/GeV) c/ VE/GeV, (3) 

where, accordingly to MC, we have fixed the constant term to 
be 5 %. We found b = 1.1% and c - 1.4% when using the 
gaussian fits to the spectra. If we repeat this procedure, for the 
fits with the logG function, we get b = 0.8% and c = 2.4%. 

We have investigated the large b/E term by measuring the 
total detector noise with a gaussian fit to Qtot in events with- 
out any electron beam impinging (see Fig. [5J. We find ctq = 
4.2H-4.8 MeV which is slightly larger than the incoherent sum 
of o~ P ed resulting to be 3.6-^3.8 MeV. A not negligible coherent 
noise is present and a much smaller (1/2) noise level has been 
previously measured in the electronic laboratory. However, the 
noise does not fully explain the large b term found. We are still 
investigating the origin of this contribution. 

We have then determined the position resolution by compar- 
ing the reconstructed centroid done with the crystals with the 
position provided by the BPM of BTF The centroids were de- 
fined as X pos - 2 QiXilQ, ot . We observe a position resolution 
of 2.8 -r 3 mm at 500 MeV, which is slightly better than the crys- 
tal pitch as expected by the application of the centroid method. 

We have finally reconstructed the calorimeter timing after 
correcting it, event by event, for the arrival time of the elec- 
trons in the LINAC spill. This was done by measuring the 
timing with the scintillators fl,f2. A very small pulse-height 
dependence on the calorimeter timing remains due to the usage 
of Constant Fraction Discriminators of KLOE. The weighted 
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Figure 7: Distribution of the average timing reconstructed with the inner matrix 
for 100 MeV electron beam. 

energy average over all calorimeter, T c ) u , was done after sub- 
tracting the average To of each cell. In Fig. [7] the distribution 
of T c i u is shown for 100 MeV electron beam. A clean gaussian 
response is observed with a time resolution, cr T , of ~ 265 ps 
at 100 MeV. A similar study at other energies provided com- 
patible results (e.g. cr T ~ 245 ps at 500 MeV). We have also 
measured the time jitter of the fl,f2 scintillators, by studying 
their time difference, which resulted to be of ~ 240 ps thus ex- 
plaining the week energy dependence of crj and pushing for a 
new measurement with a much reduced time-jitter. 
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